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I.  Problem  Definition. 

/ 

'  t  i  . 

Consider-the  problem  of  control  of  a  beam  as  shown  in  Fig.  1.  The  beam  is  moving  in  the 
x-y  plane.  It  extends  from  x=0  to  x=L.  The  left  end  at  x=0  is  clamped  to  an  actuator  which 
moves  the  beam  along  the  y-axis.  The  control  input  is  the  force  u(t)  in  y  direction.  While 
moving,  the  beam  may  vibrate.  Let  z(t)  denote  the  displacement  of  the  left  from  y=0  and 
w(t,x)  the  displacement  of  the  beam  from  the  line  y=z(t)  at  position  x  and  time  t.  Suppose  a 
position  sensor  is  place  on  the  beam  and  the  sensing  output  is  v(t,  xo)=z(t)  +  w(t,xo),  where 

i  c  ^ . 

0<xo<L  is  the  sensor  location.  We  are  interested  in  the  case  when  the  flexure  w(t,x)  of  the  beam 
is  significant.  The  problem  is  to  synthesize  a  feedback  control  taw  which  moves  the  beam  from 
one  position  to  another  in  a  stable  manner. 

It  is  well  known  [1]  that  when  the  sensor  and  the  actuator  are  colocated  a  simple  lead 
compensator  suffices  to  produce  a  stable  design.  This  result  holds  even  when  the  beam  dynamics 
are  consided  as  a  system  with  infinite  zero-damping  modes,  and  can  be  shown  using  root  locus 
arguement  [2].  This  stabilization  method  may  break  down,  however,  when  there  is  a  positional 
gap  between  the  sensor  and  actuator.  In  this  case  the  classical  compensation  techniques  are  no 
longer  effective.  Time-domain  optimization  approaches  based  on  state-space  models  have  been 

s 

applied  to  this  problem,  see  for  example  {3J.  In  this  article  vwe  present  a  case  study  of 
noncolocated  beam  control  problem  using  frequency-domain  optimization  method  proposed  by 
Professor  Kwakernaak  [4,  5]:  We  emphasize  the  choice  of  the  weighting  functions  in  the  cost 
function,  and  the  search  method  which  always  leads  to  stable  designs. 

II.  Dynamic  Models. 

Assume  the  beam  has  a  constant  mass  desity  p,  and  a  constant  bending  rigidity  El.  The 
dynamics  of  the  beam,  under  the  assumption  of  negligible  shearing,  satisfy  the  equations 


p  L  (d2  z/dt2)  +  p  02w/3t2)  dx  -  u(t) 
p  (d2z/dt2)  +  p  (<^2 w/Jt2)  +■  El  (cMw /dx4)  =  0 


(D 


with  boundary  conditions 


w(t,  0)=wx(t,  0)  =  wxx(t,  L)=  Wx x x (ti  L)  =0 


Let  <j>i  be  the  eigenmode  shapes  of  the  beam.  We  discretize  the  dynamics  by  the  decomposition 


w(t,  x)  =  £  qi(t)  <)>i(x) 


Let  £,  =  (z,  qi ,  q2,  ...  )T.  We  obtain  an  equivalent  model  of  infinite  dimension  as 


M$+K$=Fu 


y  =  c  £ 


where 


p  Widx 


P  \<t>i  «t>j  dx 


[lii± 
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F  =  (1  00  ...  )T,  and  C  =  (1  4>i(x0) <|>2(xo)  ...  ). 

Suppose  the  sensor  is  colocated  with  the  actuator  at  xo,  then  the  output  matrix  becomes 
C«(1  0  0  •••  ). 

This  input-output  relation  gives  rise  to  a  well-known  pattern  (see  [1,  61)  of  poles  and  zeros  as 
shown  in  Fig.  2.  A  simple  proof  is  given  on  the  basis  of  root  locus  arguement  in  [2]. 

Notice  that  in  this  modeling  the  beam  is  assumed  to  have  no  structural  damping,  which 
results  in  all  the  poles  and  zeros  lying  on  the  imginary  axis.  In  reality  this  is  of  course  not  the 
case.  Typically  flexibility  modes  of  very  high  frequency  tend  to  have  good  damping.  For 


Kji 


stabilizing  the  system  described  above,  it  is  therefore  possible  to  first  truncate  the  high- 
frequency  modes  and  concentrate  on  a  finite  number  of  modes  with  zero  damping. 

Consider  a  finite-dimensional  model  after  truncation.  Suppose  the  position  sensor  moves 
away  from  the  location  of  the  actuator.  It  can  be  shown  that  the  zeros  will  move  upward  along  the 
imaginary  axis.  Beginning  from  the  highest  mode  the  zeros  will  cross  the  poles  while  the  poles 
remain  at  the  same  locations.  At  some  point  the  zeros  will  break  into  complex  patterns  and  the 
system  has  nonminimum  phase  zeros.  This  observation  is  very  important,  since  nonminimum 
phase  zeros  impose  severe  constraints  to  control  synthesis.  In  Fig.  3.  the  zero  movement  due  to 
various  sensor  placements  is  shown  with  a  three  meter  auminum  beam. 

III.  Lead  Compensation. 

It  has  been  observed  (see  (1])  that  a  lead  compensator  suffices  to  stabilize  every  mode  of 
the  system  in  (4)  when  the  sensor  and  the  actuator  are  colocated.  That  is,  the  following  theorem 
holds. 

THEOREM  1 .  Consider  an  open-loop  system 

G(s)  =  Ki  n  (s  -  Zi)  /  (s2  n  (  s  -  pi)]  (6) 

Isl  It  I 

where  Ki  >  0,  and  0  <  zi  <  pi  <  zi+i  for  all  i.  Let  H(s)  be  a  lead  compensator 


H(s)  =  K2  (s/z  +  1)/(s/p  +  1),  z<p. 


(7) 


! 


Then  there  exists  a  positive  real  number  K*  such  that  for  any  K  in  the  interval  (0,  K*)  the 
closed-lop  system 


M(s)  =  G  /  (1  +GH) 


l 

(8)  • 


is  asymptotically  stable. 


A  sketch  of  the  proof  of  this  theorem  is  illustrated  in  Fig.  4.  The  departure  angle  9i  of 
the  root  locus  from  the  i-th  mode  satisfies  the  equation 

9i  =  0p  -  9z  +  ft  /  2  (9) 

where  0p  and  02  are  the  phases  of  the  pole  and  the  zero  of  the  lead  compensator.  It  is  clear  that 
for  a  small  gain  the  system  is  stable.  Now  suppose  for  a  certain  mode  the  order  of  the  pole  and 
the  zero  is  reversed  as  shown  in  Fig.  5,  which  corresponds  to  a  noncolocated  system.  With  a  lead 
compensator  the  departure  angle  from  this  pole  can  be  computed  as 

(10) 

It  is  clear  that  the  system  is  unstable  for  small  gains. 

IV.  Kwakernaak's  Minimax  Approach. 

As  shown  in  Section  III,  the  classical  compensation  techniques  are  not  effective  in  control 
synthesis  for  noncolocated  flexible  beam  with  zero-damping  modes.  Our  goal  in  this  section  is  to 
apply  the  minimax  approach  proposed  by  Professor  Kwakernaak  [4J  to  this  problem.  We 
present  only  the  procedural  aspects  of  the  design.  For  details  of  the  supporting  theorems  the 
reader  is  referred  to  the  paper  [4]. 

Consider  again  the  truncated  system  (6)  of  the  flexible  beam.  Let  y  and  <t>  be 
respectively  the  zero  polynomial  and  the  pole  polynomial.  Consider  a  unity  feedback 
configuration  and  let  H(s)  =  ^(s)/ct(s)  be  a  controller  where  C  and  a  are  polynomials.  The 
sensitivity  function  S  and  its  complement  T  are  computed  as 

S(s)  -  (1  ♦  G  H)-1  -♦o/x:  T(s)  -  1  -  S  -  vt/X  (ii) 

where  %  -  $o  +  y  C 's  closed-loop  characteristic  polymonial.  The  objective  of  Kwakernaak's 
minimax  approach  is  to  find  a  controller  H(s)  which  minimizes  the  cost  function 


Supo)  (  I  V(jca)  S(joj)|2  +  |W(ja>)  T(ju)|2  ) 


(12) 


There  V(s)  and  W(s)  are  weighting  functions.  We  denote 

Z(s)  =  V(s)  V'(s)  S(s)  S'(S)  +  W(s)  W*(s)  T(s)  T*(s)  (13) 

where  V‘(s)  =  V(-s).  The  objective  is  then  to  minimize  Supw  Z(jw).  It  is  shown  in  [4)  that  if 
an  optimizing  solution  exists,  it  is  unique  and  it  satisfies 

Z(s)  =  k2  (14) 

For  the  beam  control  problem  we  let  the  weighting  functions  be 


V(s)  =  at  /  Pi  =  (1  +  ti  s2  )  /  (t2  0  ) 

W(s)  =  -a2  /  P2  =  (  t2  Sn+2  )  /  V  (15) 

where  it,  T2  are  positive  constants  to  be  chosen  as  design  parameters.  These  weightings  satisfy 
the  conditions: 

(1)  deg(at  )  <  deg(o),  deg{a2)  =  deg(y)  +  pole  excess. 

(2)  ai,  a2.  and  Pi  have  their  roots  in  the  closed  LHP.  The  polynomials  at  and  pi  have  no 

common  roots. 

(3)  The  polynomial  y=  aiai*P;P2*  +•  a2a2*PiPt*  has  no  rootson  the  jw-axis. 

Under  these  conditions  it  can  be  shown  that  the  minimax  problem  has  a  solution.  With  these 
weightings  and  H(s)  =  (J  a,  the  condition  (14)  can  be  expressed  as 

at  ai*o  o*  *-  a2  a2*  C  C  ^  X  X*  C|6) 


Let 


— i » -#«  ■  >«i.n 


«-»  *-« 


i*%  *-<  j 


rj  =  ai  ai  *  aj  a 2* 


Consider  the  polynomial  y  -  t yX2  with  X  as  a  parameter.  Let  XQ  be  the  first  value  for  which  the 


polynomial  loses  degree  or  assumes  a  root  on  the  imaginary  axis.  It  can  be  computed  that 


Xq  =  (  Tl2  T24  )/  (tl2  +  Tl4  T24  ). 

Now  let  H(s)  be  a  reduced-order  compensator  for  which 


c  (S)  =  S  Ci  s  n-i  ,  a  (s)  =  s"+  Z  oi  s^-i 

Lsi 


With  a  particular  X  the  equation  (16)  gives  a  solution  (  £x,ox  ).  However,  a  solution  as  such 


may  be  unstable.  It  is  shown  in  [4]  that  there  exists  a  X  >  Xo,  which  yields  an  optimal 


stablilizing  compensator.  To  facilitate  searching  for  stable  solution  we  take  the  following 
approach.  Notice  that  the  equation  defining  the  closed-loop  characteristic  polynomial 


CT0  +  ^\y  =  %=  s^mf  +  Z  di  s^n-i-l 


can  be  put  into  the  form 


S(<|>,  y)  C  =  D 


where 


1  0 

a.  1 


S(o.  V) 


^  n  -  I  ^  n  -  2 


^  n  ^r*  -  t 

0  an 


0  0 


0  6,  0 

0  b2  b, 


a\  bn  bn-  ] 

«2  0  bn 

a3  0  0 


a-  0  0 


^  o 


H 


1^  l.t.f.t  *-M  4 


I .«  *-< 


with  y  =  I  bi  sn-',  6  =  sn  +  I  ai  sn‘',  C  =  (cn  •  •  •  an-i  C t  •  •  •  Cn)^".  and 

L-\ 

D  =  (dt-at  -  .  -  dn-an  dn+i  •  •  •  d2n-i)T.  The  matrix  Sfo.  y)  is  nonsigular  when  o,  and  y  are 
coprime,  and  we  have 

C  =  S-1(4>,  y)  D  (21) 

The  closed-loop  characteristic  polynomial  can  further  be  parametrized  by  a  set  of  stable  poles 


X  =  (s  +  r)  n  [  (s  +  gi)2  +  hj2  ],  r,gj  >  0.  (22) 

L  z  1 

Therefore  the  vector  D  in  (21)  can  be  replaced  by  parameters  r,  gi,  and  hi.  Using  this 
parametrization  the  search  for  the  solution  of  the  equation  ( 1 6)  is  always  pertormed  among 
stabilizing  compensators. 

V.  Numerical  Results  and  Conclusion.  . 

Consider  a  beam  control  system  with  the  rigid-body  mode  plus  a  vibrational  mode.  Its 
transfer  function  G(s)  has  poles  and  zeros  as  shown  in  Fig.  6.  Notice  the  reverse  order  of  the 

pole  and  zero.  Results  of  synthesis  with  respect  to  two  sets  of  (n,  12 )  are  shown  in  Table  1. 

With  the  stabilizing  compensator  parametrization  the  solution  of  (16)  converges  very  well.  We 
have  compared  the  controllers  obtained  minimax  approach  with  state-space  optimal  controllers. 
It  is  our  experience  that  generally  controllers  obtained  from  the  minimax  approach  have  better 
robustness. 
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Figure  1.  Beam  control  system. 
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Figure  2.  Pole-zero  pattern  of 
the  colocated  system 
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Figure  3.  Zero  pattern  of  the 
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Figure  4.  Departure  angles  of  the  Figure 

colocated  system  with 
lead  compensation. 
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Departure  angles  of  a  lead 
compensated  system  with 
poles  and  zeros  reversed. 


Figure  6.  Design  example. 


r  1=0.5  Z2  =0.07 


ri  =i-5 


Z2  =0.14 


A.  =0.1914 

closed  loop 

tf  1=21.417 

poles 

6  2=41 6.554 

6  3=-927.023 

-  2 

$  1  =-258.72 

-  2 

$  2=964.952 

-5.41 

$3=354.773 

-4.3+j3.3 

$4=223.527 

-1 .65+j5.2 

A  =0.034 

closed  loop 

4 1=24.083 

poles 

£2=552.07 
£  3=-379.09 

-0.67 

$  1  =-353.28 

-0.67 

$2=528.18 

-7.54 

$3=201.9 

-5.57+j3.24 

$4=57.44 

-2.04+j5.62 

Table  1.  Data  on  Cl, £2,  H(s)  and  closed  loop  poles 
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